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Abstract: 

We  describe  an  integrated  system  for  NC  machining  of  parts  made  of  multi-patch  surfaces. 
The  system  enables  the  user  to  generate  tool  paths  based  on  practical  metrics  such  as 
length,  curvature  and  number  of  tool  paths.  We  also  describe  two  methods  for  interfacing 
with  commercial  CAD/CAM  packages. 
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1.0  Introduction: 


Milling  is  a  commonly  used  machining  method  for  manufacturing  complex  shapes  such  as 
the  ones  found  in  the  automobile,  aerospace,  ship  building  and  health  care  industries. 
While  rapid  prototyping  is  preferred  for  manufacturing  soft  tooling,  where  the  production 
batch  sizes  are  sm^l,  numerical  control  (NC)  milling  is  the  preferred  method  to  manufac¬ 
ture  hard  tooling  which  can  be  reused  and  is  more  durable  compared  to  the  soft  tooling. 

Milling  is  the  process  of  material  removal  due  to  the  interaction  of  a  cutting  tool  moving 
relative  to  the  workpiece.  The  manufactured  part  is  obtained  as  the  tool  comes  in  contact 
with  the  nominal  designed  part  (henceforth  referred  to  as  part)  along  specific  trajectories 
called  tool  paths.  The  accuracy  and  cost  (proportional  to  the  time  of  manufacture)  of  die 
manufactured  part  are  affected  by  the  geometry  and  spacing  of  tool  paths  on  the  part 
Especially,  the  tool  paths  for  finish  machining  are  very  important,  since  finishing  cuts 
account  for  50%  of  the  total  machining  time  and  directly  affect  the  accuracy  or  surface  fin¬ 
ish  of  the  manufactured  part 

In  this  paper,  we  introduce  an  integrated  system  for  NC  machining  of  parts  made  of  multi¬ 
patch  surfaces.  The  system  enables  the  user  to  generate  tool  paths  based  on  practical  met¬ 
rics  such  as  length  and  curvature  of  the  tool  paths  and  accuracy  of  the  manufactured  part 
The  details  of  the  tool  patii  generation  scheme,  for  a  single  parametric  patch,  has  been 
described  in  [Sarma  and  Dutta,  1996a].  In  this  paper  we  focus  on  extending  the  method  to 
handle  multiple  patches,  acommon  occurence  in  CAD  models  of  industrial  parts.  For  this 
paper  we  assume  that  the  individual  patches  comprising  the  part  are  four  sided. 

This  paper  is  structured  as  follows.  The  rest  of  this  section  will  be  devoted  to  discussing 
the  various  techniques  for  tool  path  generation  being  used  in  research  and  commercial  sys¬ 
tems.  In  Section  2,  we  give  an  overview  of  our  method  for  tool  path  generation  for  a  single 
patch.  The  focus  of  this  paper  is  on  generating  tool  paths  for  parts  made  of  multi-patch 
surfaces.  This  is  described  in  Sections  3-5.  In  addition,  we  will  discuss  issues  involved  in 
the  integration  of  our  method  with  commercial  CAD/CAM  systems  in  Section  6.  We  con¬ 
clude  with  Section  7. 

1.1  Existing  Tool  Path  Generation  Techniques: 

Tool  path  generation  is  the  method  by  which  the  trajectory  of  the  tool  is  constrained  to  lie 
on  the  part.  This  has  been  traditionally  enforced  by  (a)  constraining  the  tool  path  with  the 
help  of  curves/surfaces  defined  in  the  object  space  of  the  part  and  (b)  constraining  the  tool 
path  by  specifying  curves  in  the  parametric  space  of  the  part.  Note  that  the  literature 
quoted  in  the  following  paragraphs  is  by  no  means  exhaustive  due  to  the  number  of  papers 
in  the  field  of  tool  path  generation. 

Surface  section  curves  [Bobrow,  1985]  [Huang  and  Oliver,  1994]  and  projection  curves 
[Unigraphics,  1992],  often  used  in  tool  path  generation  fall  in  the  former  category.  In  the 
latter  category  ate  iso-parametric  [Loney  and  Ozsoy,  1987]  and  iso-curvature  curves 
[Jensen  and  ^derson,  1992].  Recently,  offset  curves  on  the  part  [Suresh  and  Yang,  1994] 
[Sarma  and  Dutta,  1996a]  have  been  proposed  as  a  method  to  generate  tool  paths. 
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A  main  disadvantage  of  object  space  methods  is  that  the  user  only  has  control  over  the 
m^imum  gallop  height,  which  is  attained  at  discrete  points  on  the  manufactured  part 
[Huang  and  Ohver,  1994]  [Sarma  and  Dutta,  1996b].  Parameter  space  methods  are  advan¬ 
tageous  sinw  the  user  can  have  a  direct  control  over  the  scallop  height  of  the  manufac¬ 
tured  part  [Suresh  and  Yang,  1994]  [Sarma  and  Dutta,  1996a]. 


figure  1.  Illustration  of  Scallop  Height^  Step-Over  etc. 


■^e  methods  for  tool  path  generation  for  compound  surfaces,  usually  involve  discretizing 
the  part  mto  a  mesh  of  simpler  elements  [Lai  and  Wang,  1994]  e.g.  triangular  elements  etc. 
Ot  the  methods  mentioned  above,  [UniGraphics,  1992]  and  [Sarma  and  Dutta,  1996b]  spe¬ 
cifically  discuss  issues  related  to  multi-patch  machining.  The  extension  of  tool  path  gener- 

auon  techniques  to  multi-patch  surfaces  is  non-trivial  and  necessaiy  for  usage  in  real 
world  problems. 


1.2  Commercial  Systems  for  Tool  Path  Generation: 

Severd  commercial  systems  [UniGraphics,  1992]  [Catia,  1994]  [ProEngineer,  1995]  pro¬ 
vide  the  capability  for  tool  path  generation.  Tool  paths  used  include  plane  section  curves, 
ISO-parametric  curves,  projection  curves  and  offset  curves.  The  method  for  calculating  the 
step-over  (illustrated  in  Figure  1)  involves  specifying  (a)  the  maximum  step-over  distance 
m  object  space,  (b)  the  maximum  scallop  height  and  (c)  the  number  of  tool  paths. 

Our  method  for  tool  path  generation  [Sarma  and  Dutta,  1996a]  provides  several  new 
options  for  the  user  to  generate  tool  paths  based  on  practical  metrics  such  as  length,  curva¬ 
ture  and  number  of  tool  paths.  Such  metrics,  which  give  the  user  valuable  information 


An  Integrated  System  for  NC  Machining  of  MulU-Patch  Surfaces 


October  31.  1996 


3 


related  to  productivity  (time  and  accuracy  of  machining),  are  not  currently  available  in 
commercial  CAD/CAM  systems.  The  integration  of  our  method  for  tool  path  generation 
with  commercial  CAD/CAM  systems  is  expected  to  supplement  the  existing  methods  for 
tool  path  generation  and  provide  more  alternatives  to  the  user  for  geometrically  compli¬ 
cated  parts. 


2.0  Overview  of  Tool  Path  Generation  for  a  Single  Patch: 

For  completeness,  we  briefly  describe  our  method  for  tool  path  generation  [Sarma  and 
Dutta,  1996a]  using  offset  curves,  for  a  single  parametric  patch.  We  assume  ball  ended 
tools  and  a  3-axis  machine.  The  two  steps  involved  are:  (a)  calculating  the  step-over  and 
(b)  sequentially  generating  the  tool  paths.  The  part  is  assumed  to  be  a  parametric  patch  S: 
r(u,v)  with  unique  normals  «(«,  v).  Curves  on  the  parametric  patch  are  represented  in 
terms  of  triples  C:  {u(t),  v(t),  h(t)),  where  i^t)  and  v(t)  represent  a  parametric  curve  on 
r(u,v)  and  h(t)  is  a  distance  function. 


2.1  Offset  Method  for  Tool  Path  Generation; 

The  first  order  approximation  of  step-over  on  the  parametric  patch  S  is  calculated  by  tak¬ 
ing  offsets  of  a  given  curve  Ci  on  S.  The  step-over  calculations  can  be  symbolically 
expressed  by  (EQ  1)  and  (EQ  2).  If  C,-:  {u(t).v(t),h(t)}  is  a  scallop  curve,  then  the  tool  path 
Ci+j:  {u’(t)y(t),0]  is  obtained  using  (EQ  1).  If  Cf  {u(t)Mt).0]  is  a  tool  path,  then  the 
scallop  curve  Cj+7:  {u’(t)y(t),k(t)}  is  obtained  using  (EQ  2).  The  functions  U/t),  V/t), 
U^t)  and  V^t)  represent  a  parametric  direction  of  the  offset  The  functions  X<^(r)  and 
Xjit)  represent  the  parametric  distance  of  the  offset 


uV)\ 

v'(0/ 


U’(Q 

\vV) 


(EQl) 


(EQ2) 


We  refer  to  h(t)  as  the  scallop  height  function,  and  it  is  specified  by  the  user.  The  selection 
of  the  type  of  scallop  height  function  (i.e.  constant,  geodesic,  normal  and  distance)  is 
based  on  practical  metrics  such  as  length,  curvature  and  number  of  tool  paths.  For  more 
details  please  refer  to  [Sarma,  1996]. 

Tool  path  generation  techniques,  in  general,  attempt  to  compromise  between  conflicting 
objectives  (reducing  time  and  increasing  accuracy).  Our  method  consists  of  the  following 
sequential  steps  based  on  (EQ  1)  and  (EQ  2):  (i)  given  a  tool  path  find  the  farthest  scallop 
curve  that  satisfies  a  given  scallop  height  function  and  (ii)  given  a  scallop  curve  find  the 
farthest  tool  path  that  lies  on  the  parametric  patch.  The  procedure  begins  by  the  user  spec- 
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ifying  a  starting  curve  on  the  parametric  patch. 

Procedure: 

1.  Calculate  the  scallop  height  function  h(t)  based  on  the  local  properties  of  the  part  S 
along  the  given  starting  curve  C^. 

2.  Using  (EQ  1)  or  (EQ  2),  compute  Cy  from  Co-  If  Cq  is  a  scallop  curve  (tool  path)  then 
Cy  is  a  tool  path  (scallop  curve). 

3.  If  Cy  lies  within  the  parametric  bounds  of  S,  then  Cy, Cy  and  proceed  to  Step  1.  If  Cy 

lies  outside  the  parametric  bounds  of  S,  then  the  tool  paths  span  S  and  tool  path  genera¬ 
tion  is  complete. 

Several  parts  have  been  machined  using  this  method  and  are  reported  in  [Sarma,  1996]. 

2.2  Extensions  to  Multi-Patch  Surfaces: 

In  general,  geometrically  complex  parts  are  not  designed  with  a  single  parametric  patch. 
Hence  there  is  a  need  to  extend  the  single  patch  technique  for  tool  path  generation  to  han¬ 
dle  parts  made  of  multi-patch  surfaces.  In  addition,  we  also  consider  integration  with 
commercial  CAD/CAM  packages  and  describe  methods  for  the  same.  In  extending  our 
method  for  tool  path  generation  to  cover  parts  with  multi-patch  surfaces,  there  are  several 
issues  that  need  to  be  addressed.  They  fall  under  the  following  categories: 

1.  inputs  for  tool  path  generation 

2.  curve  data  structure 

3.  calculation  of  adjacent  curve. 

Firstly,  in  going  from  single  patch  to  multi-patch  surfaces,  the  geometry  and  topology  of 
the  part  needs  to  be  known.  This  is  address^  in  Section  3.  Next,  the  connectivity  of  the 
tool  paths  and  scallop  curves  needs  to  be  maintained  for  multi-patch  surfaces,  for  which  a 
data  structure  is  proposed.  The  details  of  the  data  structure  ate  described  in  Section  4. 
Finally,  sub-problems  related  to  calculation  of  the  adjacent  curve,  relevant  to  multi-patch 
surfaces,  need  to  be  solved.  These  sub-problems  are:  (a)  calculations  of  offsets,  (b)  calcu¬ 
lating  the  offset  across  patch  boundaries,  (c)  interpolation  of  parameter  space  curves,  (d) 
extrapolation  and  trimming  of  curves  when  necessary,  (e)  accuracy  of  object  space  curves 
and  (f)  strategy  to  deal  with  faces  that  are  not  curvature  continuous.  We  describe  tech¬ 
niques  for  these  in  Section  5. 

3.0  Inputs  for  Tool  Path  Generation: 

The  inputs  required  for  our  method  of  tool  path  generation  are:  (a)  part,  (b)  tool  radius,  (c) 
maximum  and  minimum  allowable  scallop  height,  (d)  initial  part  orientation  (important 
for  3-axis  machining),  (e)  type  of  scallop  height  function  and  (f)  starting  curve. 

The  part  is  specified  as  a  solid  or  sheet  body  [UniGraphics,  1992]  [Acis,  1995].  Bodies  are 
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the  highest  level  entities  in  geometric  modelers  which  are  able  to  keep  record  of  the  topol¬ 
ogy  and  geometry  of  the  part.  The  conceptual  structures  that  we  have  been  dealing  with 
are  as  shown  in  Hgure  2  where  we  are  interested  in  the  geometry  of  the  faces  and  edges  in 
order  to  generate  tool  paths.  We  make  the  assumption  that  all  the  faces  in  each  body  are 
joined  with  curvature  continuity  in  the  object  space.  We  discuss  exceptions  to  this  assump¬ 
tion  in  Section  5.7. 


Aas 

Body 
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1 

UNIGRAPHCIS 

1 

Shell 

Group 

1 
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Edge - 
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1 

Geometry  — Vertex 

1 

Vertex - 

Geometry 

CoEdge 

Edge - Geometry 


Vertex - Geometry 


FIGURE  2.  niustration  of  ACIS  and  UNIGRAPHICS  Data  Structures  of  Interest 

The  tool,  as  mentioned  earlier,  is  assumed  to  be  ball  ended.  The  highest  curvature  on  the 
part  corresponds  to  the  smallest  tool  radius.  A  curvature  analysis  [Maekawa  and 
Patrikalakis,  1994]  of  the  part  is  performed  to  determine  the  tool  radius.  The  maximum 
and  minimum  allowable  values  of  the  scallop  height  are  specified  by  the  user  depending 
on  the  required  accuracy  of  the  manufactured  part. 

The  initial  part  orientation  is  determined  by  examining  the  part  normals.  For  3-axis 
machining,  it  is  required  that  all  the  normals  of  the  machinable  side  of  the  part  lie  in  one 

hemisphere.  That  part  normal  which  makes  an  angle  of  no  more  than  90®  with  all  the  other 
norm^s  is  picked  as  the  tool  axis  vector.  This  is  illustrated  in  Figure  3. 


An  Integrated  System  for  NC  Machining  of  Multi* Patch  Surfaces 


October  31,  1996 


6 


FIGURE  3.  Selection  of  Tool  Axis  Given  a  Collection  of  Part  Normals 

The  type  of  scallop  height  function  (constant,  geodesic,  normal  and  distance)  can  be 
selected  depen^g  on  whether  length,  curvature  or  number  of  the  tool  paths  is  a  priority. 
The  scallop  height  function  is  switched  depending  on  the  following  criteria:  (a)  the  con¬ 
stant  scallop  height  function  is  the  default  scallop  height  function,  (b)  if  the  tool  paths 
develop  cusps  and  self-intersections  (i.e.  if  (EQ  3),  where  k  and  are  the  geodesic  cur- 

vature  of  the  tool  path  and  tool  radius  respectively,  is  close  to  being  violated)  the  geodesic 
or  normal  scallop  height  functions  can  be  used  to  compensate  and  (c)  if  the  part  geometry 
results  in  excessive  tool  path  segmentation,  the  distance  scallop  height  function  can  be 
used. 


{EQ3) 


FIGURE  4.  Selection  of  Start  Curve  for  Different  Parts 
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The  starting  curve  can  be  specified  by  object  space  constraints  such  as  (a)  constraining  the 
starting  curve  to  lie  on  the  intersection  curve  of  a  surface  and  the  part  and  (b)  constraining 
the  starting  curve  to  lie  on  the  projection  of  a  curve  (or  projection  along  a  given  direction) 
on  the  part  The  starting  curve  can  also  be  specified  by  parameter  space  constraints  such  as 
(a)  specifying  iso-parameter  curves  and  (b)  specifying  curves  in  the  parameter  space  via 
control  points.  In  our  experience,  a  boundary  curve  or  a  curve  that  is  approximately  in  the 
center  of  the  part  is  a  good  candidate  for  the  starting  curve.  The  starting  curve  could  also 
be  specified  to  be  a  closed  or  open  curve  depending  on  the  shape  of  the  part  as  shown  in 
Figure  4. 


4.0  Curve  Data  Structure: 

An  important  component  of  our  method  for  tool  path  generation  for  multiple  patches,  is 
maintaining  the  connectivity  of  the  tool  paths  and  scallop  curves  between  adjacent 
patches.  Figure  5  shows  a  part  comprising  of  a  grid  of  patches  in  the  object  space  and  in 
the  parameter  space.  For  example  the  curve  shown  in  a  dotted  line  in  Figure  5  spans  five 
surfaces.  In  order  to  traverse  the  curve,  the  “next”  and  “previous”  segments  of  Ae  curve 
have  to  be  known.  In  addition,  since  we  deal  wiA  parametric  space  curves  as  well  as 
object  space  curves,  Ae  parent  parametric  patch  has  to  be  known. 

This  mformation  is  stored  by  means  of  a  data  structure  Aat  is  schematically  shown  in  Fig¬ 
ure  6.  The  data  structure  keeps  track  of  Ae  previous  curve  (CURVE),  next  curve 
(CURVE),  parametric  space  curve  (PCURVE),  object  space  curve  (OCURVE)  and  Ae 
face  (FACE).  This  data  structure  enables  us  to  (a)  maintain  connectivity  of  tool  paths  and 
scallop  curves,  (b)  keep  information  of  parametric  curves  and  (c)  maintain  correspondence 
between  Ae  object  curves,  parametric  curves  and  faces. 

5.0  Calculation  of  Adjacent  Curve: 

In  Ais  section  we  discuss  details  of  (a)  exact  calculations  of  offsets,  (b)  calculating  Ae  off¬ 
set  across  patch  boundaries,  (c)  mterpolation  of  parameter  space  curves,  (d)  extrapolation 
and  trimming  of  curves  when  necessary,  (e)  accuracy  of  object  space  curves  and  (f)  strat¬ 
egy  to  deal  wiA  faces  Aat  are  not  curvature  continuous. 

Of  relevance  in  Ae  offset  calculation,  is  Ae  shape  of  Ae  tool  m  Ae  vicinity  of  Ae  part. 
First,  Ae  cutting  tool  is  replaced  by  a  sphere  of  Ae  same  radius.  This  simplifies  Ae  com¬ 
putation  of  Ae  swept  section  of  Ae  tool  necessary  in  our  tool  paA  generation.  The  swept 
sections  (i.e.  cross  sections  of  Ae  swept  volume)  are  circles  Aat  lie  in  Ae  normal  plane 
(i.e.  plane  spanned  by  Ae  tangent  and  normal  of  Ae  tool  paA  at  any  given  pomt)  of  Ae 
part. 


5.1  Calculation  of  Tool  Path  From  Scallop  Curve: 

Consider  a  tool  paA  C,-:  [u(t),  v(t),  0}  on  a  single  parametric  patch  S:  r(u,v).  The  goal  is  to 
calculate  Ae  scallop  curve  C/+i-  {  u'(t),  v’(t),  h(t)],  where  Ae  scallop  height  function  h(t) 
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is  known.  Let  us  consider  a  point  F,  on  C/  evaluated  at  the  parameter  tQ.  Let  a  swept  sec¬ 
tion  n  be  placed  at  P Let  the  corresponding  point  on  the  scallop  curve  be  Pg.  The  points 
P^t  ^  points  at  offsets  of  the  tool  radius  (/?)  and  scallop  height  (ho)  respectively 

from  the  parametric  patch  S.  This  is  illustrated  in  Figure  7.  Note  that  while  Figure  7  is  an 
illustration  in  two  dimensions,  the  point  Pg  does  not  necessarily  lie  in  the  plane  of  fl .  By 

definition  of  the  swept  section,  the  unit  normal  /q  of  the  plane  11  can  be  calculated  as  fol¬ 
lows: 


^uo^io  ^vO^rO 
^uO^tO  ^vO^/oj 


(EQ4) 
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FIGURE  6.  Data  Structure  for  Scallop  Curve  and  Tool  Path 

where  it  is  assumed  that  the  subscript  0  implies  evaluation  at  the  parameter  tQ  and  the  other 
subscripts  imply  partial  differentiation  with  respect  to  the  variable.  The  plane  of  the  swept 

section  is  spanned  by  the  unit  vectors  Hq  and  Bq  ,  where  Hq  is  the  unit  normal  of  the  para¬ 
metric  patch  evaluated  at  tQ  and  bo  is  a  vector  orthogonal  to  both  Hq  and  Iq  .  The  point 
can  be  found  by  evaluating  the  point  at  which  the  residual  91  shown  below  vanishes. 

91  =  /i- 

|rQ  RHq  +  RCosQbo  -  /?sin0nQ  -  ProjirQ  +  R/Iq  +  RCosQBq  -  /?sin6no)| 

The  operator  Proj{P)  represents  normal  projection  of  the  point  P  on  the  surface  5.  (EQ  5) 
can  be  solved  by  using  any  iterative  technique  e.g.  modified  regula  falsi,  etc.,  to  obtain  the 
point  Pg  and  its  parameter  values  on  the  surface  S  (the  parameter  values  are  a  by-product 
of  the  normal  projection).  The  starting  point  for  the  iterative  calculation  can  be  specified 
as  the  first  order  approximation  shown  in  Section  2.1.  It  is  possible  that  F,  and  Pg  do  not 
lie  on  the  same  parametric  patch.  Then,  the  normal  projection  of  F  is  considered  with 
respect  to  the  parametric  patch  of  Pg.  The  problem  of  detecting  the  jump  in  parametric 
patches  is  addressed  in  Section  5.3. 
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FIGURE  7.  Simplified  Illustration  of  Swept  Section  Located  on  a  Parametric  Patch 

5.2  Calculation  of  Scallop  Curve  From  Tool  Path: 

Consider  a  scallop  curve  C,-:  {u(t),  v(t),  h(t))  on  a  single  parametric  patch  S:  riu,v),  where 
the  function  h(t)  is  known.  The  goal  is  to  calculate  the  tool  path  C,+;:  {u’(t),  v’(t),  0).  Let 
us  consider  a  point  F,  on  Ci  evaluated  at  the  parameter  tg.  Let  the  corresponding  point  on 
the  tool  path  be  Pf.  The  points  P\  and  P\  are  points  at  offsets  of  the  tool  radius  (R)  and 
scallop  height  {hg)  respectively  from  the  parametric  patch  S.  Consider  a  circle  11  which  is 
placed  at  F’,  such  that  its  plane  normal  is  the  tangent  of  the  scallop  curve.  This  is  illus¬ 
trated  in  Figure  8  where,  by  definition  of  the  swept  section,  P*f  is  constrained  to  lie  on  IT 
(the  dashed  circle).  Note  that  Figure  8  is  a  two  dimensional  illustration  and  F^  does  not 
necessarily  lie  in  the  plane  of  11 . 

A  r^gUfQ  +  r^oVjo  +  +  ^o^vo^to 

*0  =  i - TT - r"7 - i  (tiQo) 

Ko“rO  +  ''voVro  +  ^0««0“r0  +  *0»vOV,o| 

_  _  ^ 

The  plane  of  the  circle  n  is  spanned  by  the  unit  vectors  hg  and  bg ,  where  Wq  the  unit 

normal  of  the  parametric  patch  evaluated  at  tg  and  Bg  is  a  vector  orthogonal  to  both  bg 

and  <0  •  The  point  F,  can  be  found  by  evaluating  the  point  at  which  the  residual  SR  shown 
below  vanishes. 
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(EQ7) 


=  /?- 

|rQ  +  /io«o  RCosQbo  +  i?sin0«o  -  PfojirQ  +  IiqHq  +  RCosQBq  +  /?sin0nQ)| 

Once  again,  (EQ  7)  can  be  solved  by  an  iterative  numerical  method  for  which  the  starting 
point  can  be  specified  by  the  first  order  approximation  given  in  Section  2.1.  It  is  possible 
that  Pf  and  Pg  do  not  lie  on  the  same  parametric  patch.  Then,  the  normal  projection  of  P  is 
considered  with  respect  to  the  parametric  patch  of  P,. 


FIGURE  8.  Simplified  Illustration  of  Swept  Section  Located  on  a  Parametric  Patch 


5.3  Offsetting  at  Patch  Boundaries: 

Figure  9  shows  a  situation  where  it  will  be  requited  to  calculate  the  offset  across  paramet¬ 
ric  patch  boundaries.  For  example,  curve  A  spans  the  parametric  patches  1, 4,  8, 9  and  6, 
while  its  offset  B  spans  1,  4,  7,  8, 9  and  6.  The  change  in  parametric  patches  is  detected 
when  the  first  order  approximation  for  calculating  offsets  (Section  2.1)  yields  an  offset 
point  outside  the  current  parametric  patch.  This  happens  when  the  parametric  coordinates 
of  the  offset  point  lie  outside  the  parametric  domain  of  the  current  patch.  There  are  three 
cases  in  which  the  offset  point  lies  outside  the  current  surface  as  illustrated  in  Figure  10. 

Case  (a)  can  be  detected  when  the  curve  A  is  locally  very  close  to  the  boundary  of  the  cur¬ 
rent  surface  S.  The  closeness  is  measured  when  parametric  distances  from  consecutive 
points  on  A  to  the  boundary  of  S  are  less  than  a  pre-specified  parametric  distance  toler¬ 
ance.  This  usually  occurs  when  the  start  curve  approaches  the  end  of  the  parametric 
domain  of  the  current  surface  S.  In  this  case,  the  edge  nearest  to  the  start  curve  A,  such  as 
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E  is  detected  and  the  surface  S’  adjacent  to  E  and  distinct  from  S  is  found.  This  surface  S 
is  then  used  as  a  candidate  for  projecting  points  on  the  offset  curve. 


FIGURE  9.  Dlustration  of  Offsets  Across  Patch  Boundaries 

Case  (b)  can  be  detected  when  the  end  points  of  curve  A  on  consecutive  surfaces  S  and  S 
are  coincident  in  the  object  domain.  This  is  the  easiest  situation  to  detect,  since  the  start 
and  end  points  of  the  start  curve  A  on  each  surface  S  and  S’  can  be  calculated  (from  the 
curve  data  structure).  In  addition,  only  one  edge  such  as  E  can  be  found  in  the  vicinity  of 
the  coincident  points  where  there  is  a  change  in  surfaces  from  S  to  S’.  In  such  a  case,  the 
surface  S’  adjacent  to  E  and  distinct  from  S,  is  used  as  a  candidate  for  projecting  the  points 
on  the  offset  curve. 

Case  (c)  can  be  detected  when  curve  A  locally  spans  three  surfaces  S,  S  and  S  .  In  addi¬ 
tion,  the  length  of  one  segment  of  curve  A  (e.g.  the  segment  of  curve  A  on  S’)  is  lesser 
than  a  pre-specified  parametric  length  tolerance.  Let  S’  be  the  current  surface.  In  such  a 
case,  the  vertex  V  closest  to  the  segment  of  curve  A  on  S’  is  found.  In  addition,  the  edges 
El  and  E2  adjoining  the  vertex  V  on  the  current  surface  S’  are  also  found.  The  candidate 
surface  for  projection  is  that  surface  which  shares  the  vertex  V  with  S’,  but  does  not  share 
edges  El  and  E2  with  S’ .  If  more  than  four  surfaces  are  connected  at  vertex  V,  then  all  the 
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surfaces  that  satisfy  the  above  condition  are  candidates  for  projecting  points. 


FIGURE  10.  niustration  of  Offsets  at  Patch  Boundaries 


5.4  Interpolation  of  Parameter  Space  Curves: 

The  result  of  offsetting  a  given  curve  on  the  parametric  patches  is  a  set  of  parameter  space 
points  with  their  corresponding  parametric  patches,  which  have  to  be  interpolated  to 
obtain  another  set  of  parametric  curves.  Since  we  are  interested  in  the  offset  property  of 
the  corresponding  object  space  curves,  we  use  a  Hermite  interpolation  of  the  parameter 
space  points  to  obtain  the  parameter  space  curves.  Hermite  interpolation  requires  the  use 
of  points  and  tangents.  Please  refer  to  any  standard  textbook  [Farin,  1990]  for  details  of 


An  Int^ratcd  System  for  NC  MactuniDg  of  Multi-Patch  Surfaces 


October  31, 1996 


14 


Hermite  interpolation.  In  this  section  we  outline  a  method  to  determine  the  tangents  at  the 
points  of  interest 

First  we  describe  the  calculation  of  tangents  of  the  scallop  curve  given  the  tangents  of  the 
tool  path.  We  refer  back  to  Figure  7  which  shows  the  swept  section  fl  with  respect  to  the 

part.  Consider  the  tangent  to  the  scallop  curve  represented  by  the  vector  To  which  will 
have  an  equation  as  shown  in  (EQ  6).  By  definition  of  the  swept  section,  the  tangents  at  the 
points  Pf  and  P\  are  required  to  be  unit  vectors  of  the  same  direction  given  by  (EQ  6). 
Hence  we  have  the  following  equation; 


f  .h  -  n  -  ( MvqSv ^  f 

®  ®  l|r„oS«+r^oSv  +  fto^„o5M+MvoH^  ° 

This  equation  can  be  solved  to  yield  the  unit  parameter  space  tangents 
the  parametric  space  tangents  have  to  be  normalized. 


(EQ8) 

as  follows,  where 


hu 

8v 


=  normalize 


-h  •  (''vo  +  Mvo) 
*o-(r„o  + V«o) 


(EQ9) 


Now  we  describe  the  calculation  of  tangents  of  the  tool  path  given  the  tangents  of  thescal- 
lop  curve.  We  refer  back  to  Figure  8  which  shows  the  circle  IT  with  respect  to  the  part 

Consider  the  tangent  to  the  scallop  curve  represented  by  the  vector  To  which  will  have  an 
equation  as  shown  in  (EQ  4).  By  definition  of  the  circle  IT ,  the  tangents  at  the  points  P\ 
and  Pj  are  required  to  be  unit  vectors  of  the  same  direction  given  by  (q  .  Hence  we  have  the 
following  equation: 


Tq  -  bo  =  0  = 


r„oSn  +  r^o6v\  . 

,r„o8w  +  r^oH^  ® 


(EQlO) 


This  equation  can  be  solved  to  yield  the  unit  parameter  space  tangents  as  follows: 


=  normalize 


f  vO 

A 

bo  • 


(EQll) 


5.5  Extrapolation  and  Trimming: 

Often  the  tool  paths  (scallop  curves)  need  to  be  trimmed  or  extrapolated  to  span  the  entire 
part  This  is  shown  in  Figure  9  where  curve  B’  has  to  be  extrapolated  at  one  end  and 
trimmed  at  the  other  end. 

The  extrapolation  is  done  by  extending  the  tangent  a  finite  amount  in  the  direction  of  the 
last  tangent  This  is  shown  in  Figure  1 1  where  an  extra  point  p  is  introduced  to  the  existing 
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set  of  points  and  tangents.  The  tangent  at  the  point  p  is  the  same  as  its  direction  of  exten¬ 
sion  i.e.  the  extrapolation  is  linear. 


\ 


FIGURE  11.  Illustration  of  Extrapolation  and  Trimmii^  in  Parameter  Space 

After  all  extrapolations  are  performed,  the  parameter  space  curve  is  obtained  by  Hermite 
interpolation.  The  resultant  curve  is  then  trimmed  by  intersecting  it  with  each  boundary  of 
the  parametric  patch.  Note  that  this  operation  is  a  curve-curve  intersection  in  a  plane.  For 
example,  in  Figure  1 1,  the  points  P  and  Q  arc  determined  by  intersecting  the  curve  inter¬ 
polated  through  the  points,  with  the  parametric  curve  boundaries. 

5.6  Accuracy  of  Object  Space  Curves: 

Deriving  the  object  space  curves  from  the  parameter  space  curves  is  necessary  since  the 
tool  paths  and  scallop  curves  are  actually  in  object  space.  In  addition,  any  subdivision  of 
the  parameter  space  curves  has  to  be  done  by  considering  the  discretization  errors  in  the 
object  space.  We  evaluate  points  on  the  object  space  curve  such  that  the  object  space  dis¬ 
cretization  error  is  lesser  than  a  maximum  allowable  value  e .  This  requires  (a)  the  calcula¬ 
tion  of  the  curvature  of  the  object  space  curve  and  (b)  evaluating  points  on  the  object  space 
curve  based  on  the  discretization  error. 

The  curvature  of  a  curve  on  a  surface  consists  of  two  components  -  the  normal  curvature 
and  the  geodesic  curvature.  The  normal  curvature  is  that  component  of  the  curvature  con¬ 
tributed  by  the  surface  on  which  the  curve  lies.  The  geodesic  curvature  is  that  component 
of  curvature  contributed  by  the  curve  meandering  in  the  tangent  space  of  the  surface  and 
hence  is  sometimes  also  referred  to  as  the  tangential  curvature. 

Consider  a  parameter  space  curve  C  represented  by  the  functions  {u(t),  v(t),  0).  Let  the 
object  space  tangent  of  the  curve  C  be  represented  by  the  vector  t .  Assume  that  the  corre- 
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spending  object  space  curve  is  parametrized  by  its  arc-length  s.  The  normal  curvature  vec¬ 
tor  K„  and  the  geodesic  curvature  vector  are  given  by  the  following  expressions 
[Struik,  1950]: 


II 

(EQ12) 

II 

X 

(EQ13) 

=  + 

(EQ14) 

where  h  is  the  unit  normal  vector  of  the  surface  and  ,  using  the  same  convention  as  in 

Section  5.1,  is  the  derivative  of  the  unit  tangent  vector  with  respect  to  the  arc  length  s.  The 
curvature  vector  of  the  curve  C  in  object  space  is  a  vector  sum  of  the  normal  and  geodesic 

curvature  vectors  as  shown  in  (EQ  14).  The  expression  for  2^ ,  in  a  form  easy  for  evaluation 
can  be  obtained  by  applying  chain  rules  of  differentiation  as  follows: 


/,  = 


2 


r,r, 


« 


(EQ15) 


FIGURE  12.  niustration  of  Discretization  Based  on  Circle  of  Curvature 
where  is  the  derivative  of  the  parameter  t  with  respect  to  the  arc  length  given  by: 


1 


h  = - 

^  r„u, +  r^v, 


(EQ16) 


The  curve  C  can  be  discretized  by  considering  the  circle  of  curvature  of  the  object  space 
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curve  as  shown  in  Figure  12  and  stepping  along  C  by  the  arc-length  a  given  by  the  follow¬ 
ing  expression: 


a 


-0  =  -acos(l-Ke) 
K  K 


(EQ17) 


5.7  C®  and  Continuous  Patches: 

In  the  case  of  parametric  patches  which  do  not  have  curvature  continuity  in  the  object 
space,  the  tool  paths  and  scallop  curves  develop  discontinuities  as  shown  in  Rgure  13. 
This  is  disadvantageous  for  two  reasons:  (a)  gouging  or  over-machining  could  occur  and 
(b)  the  length  of  the  tool  paths  could  increase  depending  on  the  connectivity  of  tool  paths. 
We  have  developed  heuristic  procedures  to  prevent  gouging  and  inefficient  tool  path  con¬ 
nectivity. 


FIGURE  13.  Illustration  of  Offsets  Across  Curvature  Discontinuous  Patches 

First  the  part  is  divided  into  regions  which  are  curvature  continuous.  If  the  curvature  con¬ 
tinuous  regions  are  large  enough  compared  to  the  part,  tool  path  generation  can  be  done 
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independently  in  each  curvature  continuous  region.  If  the  curvature  continuous  regions  are 
small  compared  to  the  part,  the  following  heuristics  are  applied. 

Case  (a):  If  adjacent  surfaces  are  curvature  continuous,  yet  yield  discontinuous  tool  paths/ 
scallop  curves  as  shown  in  Figure  14a  we  conclude  that  there  are  numerical  errors  in 
calculations.  In  such  cases  we  join  the  segments  of  the  offset  curve  as  shown  in  Figure  14a 
(right).  This  joining  is  performed  by  (i)  extrapolating  each  curve  segment  in  the  surface 
domain  such  that  it  locally  follows  the  adjacent  curve  segment  and  (iii)  trimming  the  para¬ 
metric  curve  segments. 


FIGURE  14.  Heuristics  for  Curvature  Discontinuous  Surfaces 

Case  (b):  If  adjacent  surfaces  ate  curvature  discontinuous  as  illustrated  in  Figure  4.30b 
(left)  the  offset  curves  of  A  on  each  surface  caimot  be  joined  by  interpolation  alone.  We 
propose  the  following  steps:  (a)  trim  the  offset  curves  and  (b)  introduce  an  intermediate 
curve  segment,  along  the  common  boundary  between  adjacent  surfaces,  joining  the  dis¬ 
continuous  segments  of  the  tool  path/scallop  curve  as  shown  in  Figure  4.30b  (right).  The 
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intermediate  curve  segments  need  to  be  tagged,  since  they  need  not  be  offset  any  further  - 
only  the  trimmed  curves  segments  of  B  are  candidates  for  further  offsetting. 

Case  (c):  If  adjacent  surfaces  are  curvature  discontinuous  as  illustrated  in  Figure  4.30c 
(left),  and  joining  the  corresponding  offset  curve  segments  by  intermediate  curves  causes 
significant  over-machining  on  the  common  boundary  between  the  surfaces,  the  following 
steps  are  proposed:  (a)  trim  the  offset  curves  and  (b)  introduce  an  intermediate  curve  seg¬ 
ment,  along  the  common  boundaiy  between  adjacent  surfaces,  joining  the  discontinuous 
segments  of  the  tool  path/scallop  curve  as  shown  in  Figure  4.30c  (right).  In  this  case  the 
intermediate  curve  segment  need  not  join  segments  of  the  same  offset  curve.  For  example, 
in  Figure  4.30c  (right)  an  intermediate  curve  is  introduced  between  a  segment  of  curve  C 
and  another  segment  of  curve  B. 


6.0  Integration  with  Native  System  and  Examples: 


FIGURE  !$•  Dlustration  of  Integration  with  Native  CAD/CAM  System 

Our  method  for  generating  tool  paths  has  been  implemented  on  the  Silicon  Graphics  Per¬ 
sonal  Iris  woikstation  in  C++.  The  ACIS  geometric  modeler  has  been  used  for  generating 
tool  paths.  The  user  interface  for  our  implementation  uses  the  Tcl/Tk  and  GL  libraries. 
Figure  15  illustrates  a  system  architecture  for  integration  with  native  CAD/CAM  systems 
such  as  UniGraphics,  Catia  and  ProEngineer.  In  our  case,  the  link  between  the  CAD  and 
CAM  modules  has  been  established  for  UniGraphics  by  (a)  using  a  module  built  on  the 
Acis  geometric  modeler  and  (b)  using  the  native  systems’  programming  languages. 

Figure  16  shows  an  example  of  an  intake  manifold  made  of  a  grid  of  3x2  parametric 
patches.  The  inatke  manifold  was  designed  in  the  UniGraphics  modeler.  A  maximum  seal- 
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lop  height  (/i/noji)  of  0.004  inches  and  a  minum  scallop  height  of  0.002  inches  were 
specified.  The  tool  paths  were  generated  using  UniGraphics’  GRIP  programming  lan¬ 
guage.  The  length  of  the  tool  paths  (/  inches),  the  maximum  geodesic  curvature  of  the  tool 
paths  (k  inches  and  the  number  of  tool  paths  (N)  for  different  types  of  tool  paths  are 
listed  in  Table  1  for  the  same  tool  radius  {R  inches).  The  CL  (cutter  location)  file  generated 
was  post-processed  using  the  UniGraphics  manufacturing  module. 


Type  of  Tool  Path 

1 

K 

N 

*max 

^min 

R 

Iso-Parametric 

195.3 

0.031 

34 

0.004 

- 

0.25 

Constant  Scallop 
Height 

139.9 

0.047 

33 

0.004 

0.004 

0.25 

Geodesic  Scallop 
Height 

165.4 

0.026 

36 

0.004 

0.002 

0.25 

Distance  Scall(^ 
Height 

150.4 

0.030 

33 

0.004 

0.002 

0.25 

TABLE  1.  Statistics  For  Different  Tool  Paths  for  the  Intake  Manifold 


The  tool  paths  generated  above  were  judged  with  respect  to  the  criteria  of  length,  maxi¬ 
mum  geodesic  curvature  and  number  of  tool  paths.  These  criteria  were  selected,  as 
explained  in  detail  in  [Sarma  and  Dutta,  1996a],  for  the  following  reasons:  (a)  the  time  of 
machining  is  proportional  to  the  length  of  the  tool  paths,  (b)  the  size  of  the  CL  file  and  the 
time  of  manufacturing  depends  on  the  maximum  curvature  of  the  tool  paths  and  (c)  the 
number  of  sharp  turns  and/or  tool  retractions  in  the  tool  path  is  proportional  to  the  number 
of  tool  paths.  The  constant  scallop  height  tool  paths  were  judged  to  be  most  suitable  for 
the  intake  manifold  due  to  a  good  compromise  in  the  length,  curvature  and  number  of  tool 
paths.  A  mold  was  manufactured  using  the  constant  scallop  height  tool  paths  on  a  FADAL 
3-axis  NC  milling  machine.  The  final  part  was  cast  using  an  epoxy  and  is  shown  in  Figure 
17. 


7.0  Summary: 

In  this  paper,  we  have  described  a  system  for  tool  path  generation  for  multi-patch  surfaces, 
since  parts  in  the  industry  are  usually  designed  with  multi-patch  surfaces.  We  described 
several  issues  in  tool  path  generation  that  arise  when  considering  parts  with  multiple 
patches/faces  such  as  (a)  offsetting  across  patches,  (b)  extrapolation,  interpolation  and 
trimming  of  parameter  space  curves,  (c)  accuracy  of  object  space  curves  and  (d)  curvature 
discontinuous  patches.  Each  of  these  were  implemented  using  (i)  the  Acis  geometric  mod¬ 
eler  and  (ii)  UniGraphics’  programming  language  GRIP.  We  also  describe  a  method  for 
integration  of  the  proposed  method  for  tool  path  generation  with  commercial  CAD/CAM 
systems.  Our  ongoing  work  involves  extending  our  techniques  to  patches  with  arbitrary 
trimming  curves.  Our  future  work  involves  extension  of  our  method  for  tool  path  genera¬ 
tion  for  flat  ended  tools. 
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FIGURE  17.  Intake  Manifold  Manufactured  Using  Constant  Scallop  Heights 
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